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ABSTRACT 

The sum-energy spectrum of electrons emitted in double-beta decay is a 
well-known diagnostic for the nature of the physics which is responsible for 
the decay. Three types of spectra are usually considered when these exper- 
iments are analysed: one each for the standard two-neutrino {(Ufav) decay, 
neutrinoless (PPou) and a Majoron-emitting {(3/3<p) decay. It has recently been 
shown that two other electron spectra can be possible for scalar-emitting 
modes, in addition to these traditional three. One of these is softer than 
the Standard-Model fifyv decay, while the other is intermediate between the 
fifov and the usual (3(3 V spectra. The models which predict these new spectra 
are generically more natural than those which predict the traditional (3(3^ 
spectrum, in that they can accomodate the constraints following from the 
steadily improving limits on (3(3q v decay without requiring the fine-tuning 
that is endemic to the usual models. This article reviews the properties of 
the physics which can produce the new kinds of electron spectra. 

1. Introduction 

Double-beta {(3(3) decay is an extremely rare process in which 
two nuclear neutrons simultaneously convert into two protons and two 
electrons. Within the Standard Electroweak Model (SM) this decay 
occurs at second order in the charged-current weak interactions, and 
is accompanied by the emission of two antineutrinos, giving rise to 
a characteristic {(3(32 V ) electron spectrum. Despite the extremelv long 
half-lives involved — typically 10 20 yr or more — heroic effortsu over 
the past ten years have been rewarded by its experimental detection. 

Because it is such a rare process, (3(3 decay experiments also fur- 
nish a unique window onto whatever new physics may replace the SM 
at energies very much higher than those that are directly accessible in 
today's accelerators. This is because the effects for these experiments 

*Invited talk presented to the Workshop on Double Beta Decay and Related Topics, 
Trento Italy, April 1995. 
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of new interactions can in some circumstances compete with those of 
run-of-the-mill SM decays. 

In order to be detectable in /3/3 experiments, new physics must 
have either or both of the following properties: 

1. It must violate a selection rule — e.g.: electron-number (L e ) con- 
servation — which is satisfied by the SM contribution; 

2. It must contain new particle states that are light enough to be 
produced in (3(3 decay. Since Q ~ 1 MeV is typical of the energy 
release in these decays, any such new particles must be much 
lighter than this scale. 

The purpose of this article is to outline the features of types of 
new physics which satisfy the second of these properties, producing 
(3(3 decays in which new light particles are emitted. The title refers 
to 'scalar-emitting' decay modes because, with the occasional excep- 
tion^ this new light particle is usually taken to be spinless. A spinless 
particle can be particularly well-motivated if it is a (possibly pseudo-) 
Nambu-Goldstone boson (NGB), since in that case its small mass is 
naturally understood. In particular, the focus is on comparatively re- 
cent workEra which shows that models of this sort generally have very 
different features — including qualitatively different experimental sig- 
natures, such as electron spectra — than are usually assumed in the 
analyses of current (3(3 experiments. 

2. The Electron Energy Spectrum 

The experimental quantity that is used to distinguish exotic de- 
cays from the ordinary SM events in (3(3 experiments is the shape of the 
decay rate as a function of the energies, E\ and e 2 , of the two emitted 
electrons. For instance, if no new light particles are produced, then L e - 
violating new physics can be identified if it produces decays which are 
'neutrinoless' {(3(3q v ) in the sense that only electrons emerge from the 
decaying nucleus. In this case the decay rate vanishes unless the sum 
of the two electron energies, e = E\ + e 2 , equals the released energy, Q. 

Decays into electrons plus light scalars, on the other hand, predict 
a continuous decay distribution throughout the entire kinematically 
allowed interval, 2m e < e < Q, that is distinguishable from both the 
SM Pfau distribution, and the neutrinoless (3[3 0v contribution at the 
endpoint, e = Q. Indeed, the prediction of the scalar-mediated ((3(3^) 
decay in models like the Gelmini-Roncadelli (GR) model of lepton- 
number breakingSI helped to motivate the original (3(3 experiments. 
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2.1. The Spectral Index 

The different spectra that are possible in the various decays can be 
characterized (except for (3(3q v ) in terms of a single integer p or 'spectral 
index', n. This is because the decay distribution quite generally can be 
written in the following form: 

= r (Q - e x - e 2 ) n \p x e x F{e x )] [p 2 s 2 F(e 2 )] , (1) 

where F{si) is the Fermi function which describes the distortion of 
the decay distribution due to the electrostatic field of the nucleus, and 
Pi = | pj |, for i = 1, 2, represent the magnitudes of the three-momenta of 
the electrons. We neglect the mass, /i, of the light scalar (or scalars) that 
are emitted in writing eq. ([I]) . Should [i not be negligible in comparison 
with Q, then the factor (Q — e x — e 2 ) n should be replaced by [(Q — s x — 
e 2 ) 2 — /i 2 ] n//2 . Of course, if there are several types of possible decay, the 
total spectrum is a sum of terms like eq. ([]]). 

A plot of the spectral shape which follows from eq. (jl]) for various 
choices for n, is given in Fig. 1. 

The spectral shape is determined by n because the normalization, 
r , does not depend, to a very good approximation, on the two elec- 
tron energies, e x and e 2 . This is because the most important quantity 
which sets the scale for contributions to T is the typical momenta, 
p N ~ 60 MeV, of the decaying neutrons. Since the sizes of the light 
particle momenta are set by the net energy release, Q ~ 1 MeV, they 
are negligible in their contribution to To, and in this approximation 
the spectral shape becomes determined purely by phase space and the 
Fermi functions. 

For example, the phase space of the two emitted neutrinos (plus 
the electrons) of (3(3 2v in the SM implies the corresponding spectral 
index is tism — 5. Using the phase space of a single scalar instead 
of two neutrinos similarly gives the result n GR = 1 for the spectrum 
predicted for (3(3^ decay by the GR model. With one early exception,! 
this same spectrum is also predicted by all of the alternative models for 
scalar-emitting decays which have been proposed ever since the original 
GR paper. 

Of course, Nature need not be limited to the two choices n = 5 
and n — 1, and in general other values for n might be expected to 
be possible signals for (3(3 experiments. It turns out that this naive 
expectation is true, and that models exist which predict bothBoH n = 3 
andO n = 7.Q Furthermore, some of these models — particularly those 

* As will become clear later, the cases n = 5,9, ... can also be obtained by combining 
the features which produce the n — 1,3 and 7 decays. 
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for which n = 3 — can produce observable signals in (5(5 experiments 
while preserving agreement with all other constraints, such as those 
arising from precision electroweak measurements on the Z resonance. 



Electron Spectrum 



n=7 n=5 n=3 n-1 




E = E1 + E2 



Figure 1 



The j3/3 spectrum as a function of the two electrons' total kinetic 
energy for various choices of the 'spectral index' n. n = 1 corresponds 
to the dotted line, n = 3 is the dashed line, n = 5 is the solid line and 
n = 7 is the dash-dotted line. All four curves have been arbitrarily 
assigned the same maximal value for purposes of comparison. 

2.2. How to Generate 1 

The remainder of this article is intended to briefly outline the 
properties of models with each of the values n — 1,3 and 7. Before doing 
so in detail, however, it is worth identifying the two general features 
which go into the prediction of the index n for any model.Q These are: 

1. Phase Space: As was noted above, the phase space that is asso- 
ciated with a decay into two electrons and a scalar in (5(5 V decay 
implies an index n = 1, as for the GR model. The phase space 
of each additional scalar that appears in the final state similarly 
increases n by 2, so that in the absence of other contributions to 
n, a two-scalar decay {(5(5 W ) should have: n = 3, a three-scalar 
decay: n — 5, and so on. 



4 



2. Nambu-Goldstone Bosons: As is well known, if a scalar is a NGB, 
then its couplings all must vanish in the limit of zero energy and 
momentum. This generally implies a suppression of the contribu- 
tion of such a scalar in any low-energy process, and in particular 
for (3(3 decay. Since this suppression implies that the emission 
amplitude for each NGB is proportional to a factor of the NGB 
momentum, every such particle in the (3(3 final state should also 

increase n by 2, in addition to the requirements of phase space. £ 

It is immediately clear how to generate scalar-emitting (3(3 decays 
for which n ^ 1. If two scalars whose emission amplitude is not deriva- 
tively suppressed are produced in (3(3 W decay, then the index for the 
decay should be n = 3. If N such scalars are emitted then n = 2N — 1. 
Alternatively, if a single scalar is emitted in a (3(3 V event, but this scalar 
has the derivatively-suppressed couplings of a NGB, then we again ex- 
pect n = 3. Similarly, if two such derivatively-suppressed scalars are 
emitted, then the index for the the corresponding decay is n = 7. These 
arguments are borne out by detailed calculations.^ 

2.3. A Puzzle with the GR Model 

A puzzle remains as to how the original GR model itself, and 
its many successors over the years, fit into the above counting scheme. 
After all, the light scalar which is emitted in (3(3 V decay in the GR 
model is a NGB: It is the NGB for the spontaneous breaking of lepton 
number. Yet even so, the spectral index which is predicted for this 
decay is n GR = 1, rather than n = 3 as the previous counting would 
have predicted. The resolution of this puzzle is instructive because it 
reveals an additional criterion for increasing the spectral index of a 
model. 

The resolution to this puzzleQ goes as follows. It is useful to think 
in terms of variables for which the NGB of the GR model is explicitly 
derivatively coupled. (These variables may be obtained from the stan- 
dard ones by performing a field-dependent lepton- number rotation.) In 
these new variables the NGB couples directly to the electron, as well as 
to the various neutrinos of the model. The graphs which then dominate 
(3(3^ decay at low energies turn out to be those for which the NGB of 
the final state is emitted from the external electron lines, rather than 
from the neutrino lines. 

But the emission of a massless boson with a vector coupling by an 
external electron line introduces a potential singularity into the ampli- 
tude at low momenta. Indeed if the emitted particle had been a photon 

t See, however, section 2.3 below for a qualification to this statement. 
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rather than a NGB, such graphs would really be infrared divergent. 
For NGB emission, however, the infrared singularity of these graphs is 
compensated by the derivative coupling of the massless scalar, giving a 
nonzero and finite result in the zero-momentum limit. 

In practice, then, in order to obtain a real suppression of (3(3^ 
or (3(3 W decay because of the NGB nature of the emitted scalar, it is 
also necessary to forbid the emission of the scalar from the external 
electrons in the decay. As is shown below, this is an automatic feature 
of many models once they are required to not produce experimentally 
unacceptable rates for /3/3q v decay in a natural way. 

3. Models: Preliminaries 

The remainder of the article is devoted to outlining the properties 
of the various kinds of scalar-emitting decays which the (3(3 experiments 
can see. The present section sketches those features which are required 
as preliminaries to model building, and the properties of some explicit 
models are briefly discussed in the section immediately following. Some 
of the more general conclusions that can be drawn from a comparison 
of these models are finally summarized in section 5. 

3.1. The Naturalness Problem 

The purpose of this section is to argue that a serious fine-tuning 
problem exists for virtually all of the models that have been proposed 
to date which predict n — 1 for the single-scalar (3(3^ decays, or which 
predict n = 3 for the two-scalar (3(3 W decays. This fine tuning arises 
from a naturalness issue which has strong implications for any model 
which purports to predict a detectably large scalar-emitting (3(3 decay 
rate. 

This naturalness issue hinges on the following two questions which 
any such model must address: 

1. Masses: The first question is: How can an elementary scalar have 
a mass that is smaller than Q ~ 1 MeV, and so be light enough 
to permit its production in (3(3 decay? Being some five orders of 
magnitude lighter than the electroweak scale, such a small scalar 
mass introduces a fine-tuning, or 'heirarchy', problem unless there 
is a mechanism which can protect it from virtual effects at the 
weak scale. 

2. (3(3q v Decay: The second question concerns the size that is pre- 
dicted for (3(3q v decay. A model which breaks L e generically also 
predicts a nonzero (3(3q v decay rate, and this rate must not be 
larger than the current experimental limit. Since this rate is of- 
ten related to the rate for the scalar-emitting (3(3 decay, which 
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is by assumption detectable, it can be difficult to suppress the 
one reaction without suppressing both. This makes the resulting 
constraint quite powerful. 

There are several ways in which the various existing models handle 
these issues. Two mechanisms have been proposed which can permit 
such a naturally small scalar mass. Although (somewhat surprisingly) 
supersymmetry can be used to do the jobp the resulting models are 
quite complicated and fairly contrivecLThe much simpler alternative is 
to simply follow the original workers,B£9 and to require that the light 
scalar be the NGB for an exact or approximate global symmetry. 

The second question becomes important if the symmetry for which 
the scalar is the NGB is electron number itself, as is the case for virtually 
all proposed models which predict n = 1 as the index for the scalar- 
emitting (3/3 decays. This includes essentially all models that have been 
considered until very recently, and is one of the main motivations for 
taking the models with n 7^ 1 as important alternatives. 

The problem arises because the same vacuum expectation value 
(v.e.v.), v, which breaks L e typically also generates a Majorana mass 
for the electron neutrino whose size is m UeUe ~ g e sv, where g e s is the 
relevant scalar-neutrino Yukawa coupling constant. Such a Majorana 
mass gives rise to PPq v decays which would have been seen if they 
exceed the current experimental limit, which is \m UeUe \ < 1 eV. But 
this limit cannot be satisfied simply by making g e g very small, since if 
the PP^ decat itself is to be observably large, then g eS cannot be made 
smaller thanllirll3 ~ 10 -4 . From this lower limit for g e g we learn that 
the L e -breaking v.e.v. must satisfy v < 10 keV. 

But this upper bound requires the appearance in the scalar po- 
tential of a scale that is more than 5 orders of magnitude smaller than 
the electroweak scale. As such, it poses precisely the same type of fine- 
tuning problem that would have occured if it were the scalar mass itself 
that was to be fine tuned. In either case we are led to a mass scale in 
the scalar potential that is at largest several hundred keV or so. 

It is tempting to argue that a heirarchy of the size of order 
v/M w ~ lOOkeV/lOOGeV ~ 1(T 6 is not so small if the (p - v cou- 
pling is really g e g ~ 10~ 4 , since in this case radiative corrections would 
be 5v ~ ((7efr/4vr) ~ 10~ 5 . Although seductive, this argument turns 
out to be wrong. The coupling, g e g, which appears in PP V decay is 
not simply a yukawa coupling, g, of the lagrangian, but is really an 
effective coupling in which a yukawa coupling is multiplied by a small 
mixing angle: g e ff ~ g sin 2 6. The mixing angle arises because the preci- 
sion data at LEP precludes the direct coupling of any new light scalar 
to the electroweak sector. Scalar-emitting PP decay must therefore be 
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induced through the mixing of an electroweak eigenstate with another 
state which couples to the light scalar. This mixing can occur in either 
the neutrino sector, giving sterile-neutrino mediated decays, or, for ex- 
ample, in the scalar sector. In either case, the mixing angles typically 
cannot be too big without running into conflict with other experiments. 
In sterile-neutrino models, for example, weak interaction bounds imply 
that sin# cannot be larger than 10% or so. 

Now comes the main point. Because of its suppression by powers 
of sin 9, g e g can only be as large as 10~ 4 if the underlying Yukawa cou- 
pling is considerably bigger. For example sin6 < 10% implies g > 10~ 2 . 
But it is g and not g c g which controls the radiative corrections to the 
scalar potential, so 5v is of order (g/47r) rather than (g e ff/47r) as was 
argued above. As a result the estimate using g e g underestimates the cor- 
rection to v by several orders of magnitude. More realistic estimates!™ 
show that these corrections can only be sufficiently small in restrictive 
corners of parameters space, for which new degrees of freedom (like 
sterile neutrinos) are quite light, and so are strongly constrained phe- 
nomenologically. 

To be sure, the /3j3 0u decay rate could be made naturally small 
if L e were conserved, and so if v were to vanish. But in this case the 
corresponding NGB disappears and the problem with the scalar mass 
must be solved in some other way. 

These considerations point to a natural way out of this dilemma. 
Both the scalar mass and the (3(3§ v decay rate can be naturally zero 
if: (a) the scalar is a NGB, and (b) L e is unbroken. It follows that the 
light scalar must then be a NGB for some symmetry other than that 
responsible for electron-number conservation. These two conditions are 
sufficient in themselves to imply a spectral index for the associated 
scalar-emitting (3(3 processes which is respectively n = 3 for (3(3^, and 
n = 7 for (3(3 W decay. This is easily seen since these decays can only pro- 
ceed if the emitted scalar itself carries nonzero electron number which, 
together with L e and electric-charge conservation, precludes the possi- 
bility of scalar emission from the external electron lines. 

3.2. T and Nuclear Form Factors 

Before turning to representative models for each type of decay, 
it is useful to pause to record expressions for the normalization of the 
various (3(3 decay rates in a way which facilitates the comparison of 
different models. These expressions are required in order to determine 
the kinds of couplings which would be necessary in order to obtain 
observable scalar-emitting (3(3 decay rates. 

It is useful to use the (3(3 decay rate in the form given in eq. (P, 
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with the normalization given by 



um) = {G llZy )A WA)|2 ' (2) 

where G F is the Fermi constant, 9 C the Cabibbo angle, and A an am- 
plitude which depends on the decay process being computed (ft Pi, with 
% = 2u, (f or <fip), on the couplings of the model, and on some soon-to- 
be-identified nuclear matrix elements. 

A can be written, for + — > + transitions, as a convolution: 

a = j^irwwM, (3) 

where depends on the detailed properties of the leptons in the 
model, and W^i) 
vant to the decay: 



model, and W^ v {i) contains the nuclear matrix elements that are rele 



{i) = (27r)3 V wi' I d ' x {N ' lT * Mx)J » m lN) eltx - (4) 



Here J M = W7 M (1 + 7s)rf is the weak charged current that causes transi- 
tions from neutrons to protons, and \N) and \N') represent the initial 
and final + nuclei in the decay. E and M are the energy and mass of 
the initial nucleus, N, while E' and M' are the corresponding properties 
for the final nucleus, N'. The symmetries^ of the problem ensure that 
the most general possible form for isa: 

+w 5 {tfiUv - tvU^) + iw 6 e fiUap u a £ p , (5) 

where w M is the four-velocity of the initial and final nucleus, and the 
six Lorentz-invariant form factors, w a = w a (u ■ £, £ 2 ), a = 1, . . . , 6, are 
functions of the two independent invariants that can be constructed 
from £ u and u u . 

These form factors can be relatedH to the nuclear matrix elements 
as they are quoted in the literature .Iliril3 For example, in many situations 
(such as PPiv, PPou and some kinds of PP V and PP W decays) U* v oc rf v 
and so the decay rate only depends on the combination W^^. This is 
simply the difference between the Fermi and Gamow- Teller form fac- 
tors, which are defined by w F = W 00 and w GT = J2l=i Wkk- These form 
factors, when computed using the closure and nonrelativistic impulse 
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approximations in a model of the nucleus, become 
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Here e = -E — M is the average excitation energy of the intermediate 
nuclear state, r nm is the separation in position between the two decaying 
nucleons, 1 is the spatial component of in the nuclear rest frame, 
and g v ~ 1 and g A ~ 1.25 are the vector and axial couplings of the 
nucleon to the weak currents. Finally, ((N'\0\N)) represents a reduced 
matrix element from which the nuclear centre-of-mass motion has been 
extracted, and r+ (or a n ) are the raising operators for nuclear isospin 
(or the nuclear spin operators) acting on the n'th nucleon. 

4. Models: Specific Examples 

We now turn to the details of representative models for which 
n = 1, 3 and 7. Although much of what follows applies more generally 
we choose, for simplicity and for the purposes of comparison, models 
in which the corresponding (3(3 decay proceeds due to the mixing of 
the electron neutrino, u e , with a collection of sterile neutrinos, iVj. The 
reader who is not interested in the specific properties of these models, 
can skip to the next, concluding, section in which their general features 
are compared. 

4-1. The Case n — 1 

The simplest models to build are those for which the spectral 
index takes its traditional value, n = 1. Theories of this sort, when 
constructed usine sterile neutrinos, are similar to the original singlet- 
majoron models.EI They, as well as other variants with n = 1 which 
do not rely on sterile neutrinos to produce (3(3 V decays,t3 have recently 
received renewed attention. 

Suppose the neutrino mass eigenstates are denoted generically by 
Vi, and their overlap with the electron- neutrino flavour eigenstate is 
called V e i. Suppose also that represents the light scalar of the model. 
Then, taking the general Yukawa coupling lagrangian between these 
particles to be 

£ W v = -- h{a>ijlL + bijlii) Vj f* + c.c. , (7) 
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and evaluating the Feynman graph of Fig. 2, gives an amplitude,B A of 
the form of eq. @, withf] 



1.1 



(aijmjmj - £ 2 bjj) rf^ 
(P + mf -ie)(£ 2 + m 2 - is) 



(8) 





Figure 2 



The Feynman graph which is responsible for (3(3^ decay in models for 
which (3(3 decay can arise because of sterile-neutrino exchange. 

The remaining question is whether an observably large f3f3 v decay 
rate can be obtained from these expressions without becoming in con- 
flict with any existing phenomenological bounds. This can be addressed 
only by constructing an explicit model in which all of the relevant ob- 
servables can be computed as a function of a common set of underlying 
parameters. We therefore next present an representative model for this 
kind of theory. 

An explicit sterile-neutrino model which produces (3(3^ decays 
with n = 1 is simple to construct B Add to the SM two electroweak- 
singlet, left-handed neutrinos s ± , which carry ±1 unit of an unbroken 
global lepton number symmetry. Add also a singlet scalar field with 
lepton number —2. The most general renormalizable couplings of these 



tOur conventions here use rf 1 " = diag( — V ++). 
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new particles, consistent with lepton number conservation and the SM 
gauge symmetries, are 

C = -\LH^ R s_ - Ms + 7 fl s_ - ~ y ^- TrS ~ V9 * + C ' C ' ( 9 ) 

Here L and H respectively denote the usual SM lepton- and Higgs 
doublets, and j L and ^ R denote the usual projections onto left- and 
right-handed spinors. Majorana spinors are used throughout to rep- 
resent the neutrinos, so the spinor conjugate used above employs the 
charge-conjugation matrix, C, according to I/j = vJC~ x . 

As discussed previously, the experimental absence of (3(3q v decay 
requires the expectation value of the scalar field, (</?), not to be larger 
than ~ 100 keV. Rather than fine-tuning in such a small scale, it is sim- 
pler to take (</?) = 0, so that lepton number is not spontaneously broken. 
Then the fine-tuning simply becomes the requirement that the scalar 
mass be < 10 keV. In this case, the spectrum contains three massless 
neutrinos, z/ e , and z/ T , together with a massive Dirac neutrino, z/^, 
which mixes with the electron-type charged-current weak interactions. 

The model can have a detectable^ (3(3^ decay rate, provided that 
the masses and mixings are chosen judiciously. Because the scale of 
the factor W^ v in eq. (|3|) is set by the nucleon momentum, p N ~ 60 
GeV, (3(3 decay rate tends to become suppressed if all neutrino masses 
are much larger than, or much smaller than this scale. As a result, 
a successful choice of parameters, which can also avoid bounds from 
other laboratory experiments, puts the heavy neutrino mass eigenstate 
at several hundred MeV, with a fairly large (somewhat less than 10%) 
mixing with v e . More of the qualitative features and problems that arise 
with these models are outlined in section 5, below. Those interested in 
the details of the analysis are referred to the recent literature .H'ta 

4-2. The Case n = 3: Two-scalar Decays 

There are two possibilities for producing n = 3 decays: two light 
scalars could be emitted without a NGB suppression of the emission 
amplitude, or one light NGB-suppressed scalar could be emitted. An 
example of the two-scalar decay is given here, even though the two- 
scalar emission models involve the same fine-tuning problems as do the 
n = 1 models just described. The alternative models, for which n = 3 
arises in single-scalar decays, are the subject of the next subsection. 

Consider in this theory of two types of left-handed neutri- 

nos, Ui and iV , which respectively carry lepton number L e {y() = +1 
and L e (N a ) = 0, and which are coupled to the light scalar, ip, which 
has lepton number L e (ip) = +1. The most general renormalizable and 
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L e -conserving Yukawa couplings involving these fields are: 



C 



yuk 



Vi (A ia -f L + B ia ^ R ) N a ip + h.c, 



(10) 



where A ia and B ia represent arbitrary Yukawa-coupling matrices. These 
neutrinos are endowed with a set of generic lepton-number-conserving 
masses, m u . and m Na . For the L e = neutrinos, N a , this is accomplished 
by simply introducing a general majorana-mass term. For the L e = 1 
neutrinos, however, an additional collection of L e = 1 right-handed 
neutrinos are required, with which L e -invariant Dirac masses may be 
formed. 

This type of theory produces decay due to the Feynman 

graph of Fig. 3. Two scalars must be emitted in this decay because 
of L e conservation. Evaluating this graph gives a result of the form of 
eqs. H and || withH: 



V eVi V ev M ja rT 



(£ 2 + ml. - ie) (I 2 + m 2 . - ie) (£ 2 + m 2 Na - ie) 



where the factor, Afij a , denotes: 



'IF 



Nija = (—£ ) Ai a Bj a m Vi -^Aj a Bi a m v .-\-Bi a Bj a m Na -\-Ai a Aj a m Ui m U:j m Na . 

(12) 



N 




<P 




Figure 3 



The Feynman graph which is responsible for {3f3 vv decay in models for 
which (3/3 decay arises because of sterile-neutrino exchange. 
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It is just possible to obtain a detectably large decay rate in this 
kind of model without running into conflict with other laboratory bounds. 
Because of the comparatively soft electron spectrum (since n = 3), the 
total integrated decay rate tends to be suppressed compared to n = 1 
models by an additional two powers of the small ratio Q/p N ~ 1CT 1 . 
The additional phase space also introduces additional suppression due 
to dimensionless factors of 1/2%. As a result the total rate tends to be 
much smaller than in a comparable model for which n = 1. This makes 
it more difficult to obtain observable decays, and a sufficiently large 
rate generally requires sterile neutrinos to lie in the mass range in the 
vicinity of 100 MeV which optimizes the (3j3 decay rate. 

4-3. The Case n = 3: Single-scalar Decays 

We next turn to models which predict only single-scalar (3(3^, but 
for which the spectrum nevertheless has n = 3. This may be ensured 
by constructing the model to ensure that the emitted scalar is a NGB 
which carries two units of conserved electron number. 

The rate for (3j3 v decay is given by evaluating the result for the 
Feynman graph of Fig. 2. using the generic Yukawa coupling of eq. (0). 
By virtue of the light scalar being a NGB carrying unbroken electron 
number, one finds that the leading result in powers of the lepton mo- 
menta, as given by expression |8|, vanishes identically. It is therefore 
necessary to work to next order in these momenta, which raises the 
resulting spectral index from n = 1 to n — 3. The resulting amplitude 
then satisfies It*" = -L uf ", with!: 

L 0m = -4V2j2V et V ej h 

ij 

L mn = _ A ^2j2V ei V eJ h 



(£ 2 + m 2 - ie){£ 2 - m 2 + is) 



(£ 2 + mf — ie)(£ 2 — m 2 + is) 



(13) 



An example of a renormalizable model for which this is the decay 
formula is givenH by requiring the theory to have a nonabelian flavour 
symmetry, G = SU F (2) x U L /(1) which gets broken down to the unbro- 
ken electron number. To implement this symmetry-breaking pattern, 
introduce an electroweak-singlet scalar field, $j, which transforms un- 
der G like (2,-1). Also introduce the electroweak-singlet left-handed 
neutrino fields, N ~ (2, 0) and s ± ~ (1, ±1). 

The most general renormalizable lagrangian involving the new 
fields which respects all of the symmetries is 

C = -\LH lR s. - Ms +lR S- - g+ (N lL s+) $ - (A 7i s_) $ + c.c. 

(14) 
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$ = ir 2 $* represents the conjugate SU F (2) doublet, where r 2 is the 
second Pauli matrix acting on flavour indices. The scalar potential is 
then chosen to ensure that $ gets a VEV, which we assume to take the 

form ($) = g). 

The soft n = 3 spectrum suppresses the integrated decay rate, 
giving a detectable spectrum only when all couplings and masses are 
optimal. This once again places the sterile neutrinos in the mass range 
of several hundred MeV, with significant mixing with the electron neu- 
trino. 

4-4- The Case n = 7 

As a final example, consider the case with spectral index n = 7. 
This decay is produced by ensuring that the light scalar is a NGB, and 
that it carries only one unit of electron number so that the decay rate is 
derivatively suppressed, and two scalars must be emitted in the decay. 

The 13 (3 w decay rate is found by evaluating the Feynman graph 
of Fig. 3. Once again, the symmetry properties make it necessary to 
work to higher order in the lepton momenta than was necessary for the 
n = 3 PPtptp decay. It turns out that the expression for the resulting 
amplitude is quite cumbersome when given in terms of the couplings 
of eq. ([U]), and so we use instead variables for which the derivative 
coupling nature of the Goldstone bosons is manifest from the outset. 
The trilinear coupling to neutrinos of a Goldstone boson carrying L e = 
1, then becomes: 

£ g b = -i Va^(X ia ^ L + Y ia rf R )N a d^(p + h.c, (15) 

where the coefficients X ia and Y^, are coupling matrices that can be 
computed in any specific model.&B Using these interactions to evaluate 
the amplitude given by the diagram of Fig. 3 gives, for the special caseO 
Yia = 0: 



l^(£) = ( -4s)' y 

K 1 vio5ttV t-r 



i j a 



V eVl V ev .N ija r) 



(£ 2 + ml. — ie) (£ 2 + m 2 . — ie) (£ 2 + m 2 N — ie) 

"(16) 

Afija represents the following expression: 

Af ija = (~£ 2 )X ia X ja m Na . (17) 

It is fairly unlikely that this decay will be found in j3/3 experiments 
in the forseeable future, even if Nature should actually work contain 
n = 7 ftfltptp decays. There are two reasons why this is so. First, the 
very high spectral index for this decay, n = 7, raises several problems 
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for detecting this kind of decay. Firstly, since the electron spectrum 
is softer than that of the SM fifiiv decay, the decay electrons tend to 
come out with low energies, and so would be difficult to distinguish 
from background. 

Secondly, the soft spectrum greatly suppresses the integrated total 
decay rate, making it very difficult to get a detectable decay in a model 
which also satisfies all other laboratory bounds. Even though models 
which produce this spectrum have been constructed,u none have been 
found which are both phenomenologically viable and which predict a 
detectable (3(3 W decay rate. 

5. Models: General Features and Conclusions 

A comparison of models in these four classes leads to a number 
of reasonably robust conclusions. 

1. The models for which n — 1 and those which predict n = 3 
(3(3 W decays illustrate in detail the general fine-tuning problem 
that was argued in section 2 to be endemic to these kinds of 
(3(3^ decay. The requirement that a light scalar should exist, and 
that (3(3q v should not be predicted at an unacceptably large rate, 
taken together require a fine tuning of the scalar potential to 
ensure either an extremely small scalar mass, or an equally small 
L e -breaking v.e.v.. Although supersymmetric models along these 
lines have been constructed^ which are technically natural, they 
are also quite contrived and complicated. 

2. Models with softer electron spectra give smaller integrated decay 
rates, since the total decay becomes suppressed by higher powers 
of the small endpoint energy Q. This implies that theories pre- 
dicting n — 1 decays have an easier time producing observably 
large decay rates than do models for which n is larger. As a re- 
sult these models tend to offer the largest latitude to accomodate 
other laboratory limits and phenomenological constraints. 

To produce acceptably large scalar-emitting (3(3 decay rates, mod- 
els for which n = 3 must typically have all of the relevant dimen- 
sionless couplings be 0(1), have the mixings of the relevant sterile 
neutrino be as large as are phenomenologically allowed (< 10%), 
and have the participating new neutrino states have masses in 
the optimal mass range of a few hundred MeV. This mass range 
is prefered since it is comparable to the typical momenta, p N of 
the decaying nucleons within the nucleus, and so does not lead to 
suppressions of the form of M/p N or p N /M. The resulting models 
therefore can work, but do not leave a great deal of freedom to 
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accomodate other limits. Of the n = 3 models, those which emit 
only a single scalar tend to predict larger decay rates than the 
two-scalar decays because they are not suppressed by additional 
small phase-space factors. 

Models for which n = 7 are the worst case, and have decay rates 
that are sufficiently suppressed by powers of Q/p N that they are 
unlikely to be experimentally detectable for the forseeable future. 

3. Models with electron spectra as soft as the n = 7 decays are 
also harder to detect for another reason, besides the size of their 
total decay rate. Most of the background in (3(3 decay experiments 
occurs for electron energies which are comparatively soft, and so 
it is the soft electrons which are the ones which are hard to dig 
out of this background. This makes it all the more unlikely that 
these decays will turn up in the experimental data. 

4. Because the contribution of sterile neutrinos to scalar-emitting 
(3(3 decays become suppressed when their masses are much larger 
or much smaller than around 100 MeV, and since couplings and 
mixings tend to have to be large in order to produce an observable 
decay, the experimental discovery of (3(3 V or /3(3 VV decay would 
strongly suggest the existence of new particles in this mass range. 
The signals for such particles could be: anomalous bumps in K — * 
ev or 7r — ► ev decays; violations of weak universality in leptonic 
it decays, possible Zenlike monojet events at LEP, etc.. Signals 
in beam dump experiments would not necessarily be expected, 
since the heavy neutrinos would dominantly decay invisibly into 
ordinary neutrinos and the light scalars. 

5. All models which produce scalar-emitting (3(3 decays necessarily 
have at least one very light scalar which is significantly coupled 
to the electron neutrino. As a result, all of these models generi- 
cally run into trouble with Big Bang nucleosynthesis. The scalars 
tend not to decouple from the ordinary neutrinos, and so tend 
to violate the constraints on the gravitating degrees of freedom 
which can exist at the epoch around T ~ 1 MeV. This bound 
can be evaded for comparatively special values of the masses and 
couplings of the particles involved.cSE§u The loophole arises since 
there can be an interval during which the neutrinos are no longer 
in chemical equilibrium with the protons and neutrons, but where 
the n/p ratio has still not frozen out. During this interval the 
annihilation of sterile neutrinos can raise the electron neutrino 
abundance, which in turn acts to deplete the neutron abundance. 
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A sterile neutrino which annihilates in this way effectively counts 
as a negative number of neutrino species, and so can counteract 
the positive contribution of the light scalars. 

Of course, the details of the loophole are less interesting than 
the simple fact that the loophole exists. Although nucleosynthesis 
considerations generically disfavour models with additional light 
scalars, nucleosynthesis should and would be re-evaluated should 
a scalar-emitting (3(3 decay mode be observed. 

6. The nuclear matrix elements that are required to evaluate the (3(3 
decay rates in essentially all of the models are the same as those 
which have long been studiedE3ii3liZl within the context of (3(3q v and 
(3(3^ decay in the GR model. This can be seen since L^ u oc rf v , 
which implies that the decay rate depends only on the nuclear 
form factor combination: W^^ = w GT — w F . This is precisely the 
same combination as appears in Pfou and (3{3o v decays, and so it 
is well studied in the literature. 

The exception to this rule are those models which predict n = 3 
single-scalar (3(3^ decay, which instead depend on the antisym- 
metric part of W^. The corresponding matrix elements are less 
well studied, with a corresponding greater uncertainty in the pre- 
dicted decay rates. (Explicit formulae for these matrix elements 
in terms of nucleon operators are given in the literature .□) 

7. All of the models given here predict the same decay distribution 
as a function of the opening angle of the two electrons, so this 
observable cannot be used to distinguish one from another (as it 
can be used to distinguish decays mediated through right-handed 
currents). The only exception to this statement among the models 
considered are those for which the final electrons are emitted from 
an L e = 2 scalar having electric charge Q = —2. Unfortunately 
the bounds on the masses of such scalars make the resulting (3(3 
decay rate undetectable. Q 

To summarize, (3(3 experiments can legitimately expect to see 
scalar emitting decays even though the original models which proposed 
this decay mode have since been ruled out by constraints such as those 
coming from LEP. If such decays are seen, they are most likely to point 
to new physics with properties which are quite different than what 
would have been expected from these traditional models. Although de- 
cay spectra with n — 1,3, 5, 7, ... are all possible, the basic combinations 
are n — 1,3 and 7. The n = 7 spectrum is very unlikely to seen, however, 
due to the strong suppression of this decay rate by powers of the small 
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energy release, Q. Two-scalar n = 3 decays also tend to be suppressed 
compared to single-scalar decays having the same spectrum due to the 
presence of additional small phase-space factors. Furthermore, natural- 
ness considerations (driven by the strong bound on the occurence of 
the flfl Qv decay mode) disfavour those models which predict n = 1 flfl v 
decays, or n = 3 flfl w decays. Thus, from a theoretical perspective the 
n — 3 (3 fly deays are the most natural. If such decays are seen we may 
also probably expect some new developments in nucleosynthesis and in 
precision experiments such as those which constrain lepton universality. 
May we live to see such exciting times! 
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